Vitamin B12 variants structure soil microbial communities despite soil's vast reservoir of B12.

Soil microbial communities perform critical ecosystem services through the collective metabolic activities of numerous individual organisms. Most microbes use corrinoids, a structurally diverse family of cofactors related to vitamin B12. Corrinoid structure influences the growth of individual microbes, yet how these growth responses scale to the community level remains unknown. Analysis of metagenome-assembled genomes suggests corrinoids are supplied to the community by members of the archaeal and bacterial phyla Thermoproteota, Actinobacteria, and Proteobacteria. Corrinoids were found largely adhered to the soil matrix in a grassland soil, at levels exceeding those required by cultured bacteria. Enrichment cultures and soil microcosms seeded with different corrinoids show distinct shifts in bacterial 16S composition, supporting the hypothesis that corrinoid structure can shape communities. Environmental context influenced both community and taxon-specific responses to specific corrinoids. These results implicate corrinoids as key determinants of soil microbiome structure and suggest that environmental micronutrient reservoirs promote community stability.

Rhizosphere and detritusphere habitats modulate expression of soil N-cycling genes during plant development.

IMPORTANCE: Plant roots modulate microbial nitrogen (N) cycling by regulating the supply of root-derived carbon and nitrogen uptake. These differences in resource availability cause distinct micro-habitats to develop: soil near living roots, decaying roots, near both, or outside the direct influence of roots. While many environmental factors and genes control the microbial processes involved in the nitrogen cycle, most research has focused on single genes and pathways, neglecting the interactive effects these pathways have on each other. The processes controlled by these pathways determine consumption and production of N by soil microorganisms. We followed the expression of N-cycling genes in four soil microhabitats over a period of active root growth for an annual grass. We found that the presence of root litter and living roots significantly altered gene expression involved in multiple nitrogen pathways, as well as tradeoffs between pathways, which ultimately regulate N availability to plants.

A subset of viruses thrives following microbial resuscitation during rewetting of a seasonally dry California grassland soil.

Viruses are abundant, ubiquitous members of soil communities that kill microbial cells, but how they respond to perturbation of soil ecosystems is essentially unknown. Here, we investigate lineage-specific virus-host dynamics in grassland soil following "wet-up", when resident microbes are both resuscitated and lysed after a prolonged dry period. Quantitative isotope tracing, time-resolved metagenomics and viromic analyses indicate that dry soil holds a diverse but low biomass reservoir of virions, of which only a subset thrives following wet-up. Viral richness decreases by 50% within 24 h post wet-up, while viral biomass increases four-fold within one week. Though recent hypotheses suggest lysogeny predominates in soil, our evidence indicates that viruses in lytic cycles dominate the response to wet-up. We estimate that viruses drive a measurable and continuous rate of cell lysis, with up to 46% of microbial death driven by viral lysis one week following wet-up. Thus, viruses contribute to turnover of soil microbial biomass and the widely reported CO2 efflux following wet-up of seasonally dry soils.

Expression of macromolecular organic nitrogen degrading enzymes identifies potential mediators of soil organic N availability to an annual grass.

Nitrogen (N) is frequently limiting to plant growth, in part because most soil N is present as polymeric organic compounds that are not readily taken up by plants. Microbial depolymerization of these large macromolecular N-substrates gradually releases available inorganic N. While many studies have researched and modeled controls on soil organic matter formation and bulk N mineralization, the ecological-spatial, temporal and phylogenetic-patterns underlying organic N degradation remain unclear. We analyzed 48 time-resolved metatranscriptomes and quantified N-depolymerization gene expression to resolve differential expression by soil habitat and time in specific taxonomic groups and gene-based guilds. We observed much higher expression of extracellular serine-type proteases than other extracellular N-degrading enzymes, with protease expression of predatory bacteria declining with time and other taxonomic patterns driven by the presence (Gammaproteobacteria) or absence (Thermoproteota) of live roots and root detritus (Deltaproteobacteria and Fungi). The primary chitinase chit1 gene was more highly expressed by eukaryotes near root detritus, suggesting predation of fungi. In some lineages, increased gene expression over time suggests increased competitiveness with rhizosphere age (Chloroflexi). Phylotypes from some genera had protease expression patterns that could benefit plant N nutrition, for example, we identified a Janthinobacterium phylotype and two Burkholderiales that depolymerize organic N near young roots and a Rhizobacter with elevated protease levels near mature roots. These taxon-resolved gene expression results provide an ecological read-out of microbial interactions and controls on N dynamics in specific soil microhabitats and could be used to target potential plant N bioaugmentation strategies.

Propagation of viral genomes by replicating ammonia-oxidising archaea during soil nitrification.

Ammonia-oxidising archaea (AOA) are a ubiquitous component of microbial communities and dominate the first stage of nitrification in some soils. While we are beginning to understand soil virus dynamics, we have no knowledge of the composition or activity of those infecting nitrifiers or their potential to influence processes. This study aimed to characterise viruses having infected autotrophic AOA in two nitrifying soils of contrasting pH by following transfer of assimilated CO2-derived 13C from host to virus via DNA stable-isotope probing and metagenomic analysis. Incorporation of 13C into low GC mol% AOA and virus genomes increased DNA buoyant density in CsCl gradients but resulted in co-migration with dominant non-enriched high GC mol% genomes, reducing sequencing depth and contig assembly. We therefore developed a hybrid approach where AOA and virus genomes were assembled from low buoyant density DNA with subsequent mapping of 13C isotopically enriched high buoyant density DNA reads to identify activity of AOA. Metagenome-assembled genomes were different between the two soils and represented a broad diversity of active populations. Sixty-four AOA-infecting viral operational taxonomic units (vOTUs) were identified with no clear relatedness to previously characterised prokaryote viruses. These vOTUs were also distinct between soils, with 42% enriched in 13C derived from hosts. The majority were predicted as capable of lysogeny and auxiliary metabolic genes included an AOA-specific multicopper oxidase suggesting infection may augment copper uptake essential for central metabolic functioning. These findings indicate virus infection of AOA may be a frequent process during nitrification with potential to influence host physiology and activity.

Quantitative Stable-Isotope Probing (qSIP) with Metagenomics Links Microbial Physiology and Activity to Soil Moisture in Mediterranean-Climate Grassland Ecosystems.

The growth and physiology of soil microorganisms, which play vital roles in biogeochemical cycling, are shaped by both current and historical soil environmental conditions. Here, we developed and applied a genome-resolved metagenomic implementation of quantitative stable isotope probing (qSIP) with an H218O labeling experiment to identify actively growing soil microorganisms and their genomic capacities. qSIP enabled measurement of taxon-specific growth because isotopic incorporation into microbial DNA requires production of new genome copies. We studied three Mediterranean grassland soils across a rainfall gradient to evaluate the hypothesis that historic precipitation levels are an important factor controlling trait selection. We used qSIP-informed genome-resolved metagenomics to resolve the active subset of soil community members and identify their characteristic ecophysiological traits. Higher year-round precipitation levels correlated with higher activity and growth rates of flagellar motile microorganisms. In addition to heavily isotopically labeled bacteria, we identified abundant isotope-labeled phages, suggesting phage-induced cell lysis likely contributed to necromass production at all three sites. Further, there was a positive correlation between phage activity and the activity of putative phage hosts. Contrary to our expectations, the capacity to decompose the diverse complex carbohydrates common in soil organic matter or oxidize methanol and carbon monoxide were broadly distributed across active and inactive bacteria in all three soils, implying that these traits are not highly selected for by historical precipitation. IMPORTANCE Soil moisture is a critical factor that strongly shapes the lifestyle of soil organisms by changing access to nutrients, controlling oxygen diffusion, and regulating the potential for mobility. We identified active microorganisms in three grassland soils with similar mineral contexts, yet different historic rainfall inputs, by adding water labeled with a stable isotope and tracking that isotope in DNA of growing microbes. By examining the genomes of active and inactive microorganisms, we identified functions that are enriched in growing organisms, and showed that different functions were selected for in different soils. Wetter soil had higher activity of motile organisms, but activity of pathways for degradation of soil organic carbon compounds, including simple carbon substrates, were comparable for all three soils. We identified many labeled, and thus active bacteriophages (viruses that infect bacteria), implying that the cells they killed contributed to soil organic matter. The activity of these bacteriophages was significantly correlated with activity of their hosts.

Substantial differences in soil viral community composition within and among four Northern California habitats.

Viruses contribute to food web dynamics and nutrient cycles in diverse ecosystems, yet the biogeographical patterns that underlie these viral dynamics are poorly understood, particularly in soil. Here, we identified trends in soil viral community composition in relation to habitat, moisture content, and physical distance. We generated 30 soil viromes from four distinct habitats (wetlands, grasslands, woodlands, and chaparral) by selectively capturing virus-sized particles prior to DNA extraction, and we recovered 3432 unique viral 'species' (dsDNA vOTUs). Viral communities differed significantly by soil moisture content, with viral richness generally higher in wet compared to dry soil habitats. However, vOTUs were rarely shared between viromes, including replicates <10 m apart, suggesting that soil viruses may not disperse well and that future soil viral community sampling strategies may need to account for extreme community differences over small spatial scales. Of the 19% of vOTUs detected in more than one virome, 93% were from the same habitat and site, suggesting greater viral community similarity in closer proximity and under similar environmental conditions. Within-habitat differences indicate that extensive sampling would be required for rigorous cross-habitat comparisons, and results highlight emerging paradigms of high viral activity in wet soils and soil viral community spatial heterogeneity.

Methane-derived carbon flows into host-virus networks at different trophic levels in soil.

The concentration of atmospheric methane (CH4) continues to increase with microbial communities controlling soil-atmosphere fluxes. While there is substantial knowledge of the diversity and function of prokaryotes regulating CH4 production and consumption, their active interactions with viruses in soil have not been identified. Metagenomic sequencing of soil microbial communities enables identification of linkages between viruses and hosts. However, this does not determine if these represent current or historical interactions nor whether a virus or host are active. In this study, we identified active interactions between individual host and virus populations in situ by following the transfer of assimilated carbon. Using DNA stable-isotope probing combined with metagenomic analyses, we characterized CH4-fueled microbial networks in acidic and neutral pH soils, specifically primary and secondary utilizers, together with the recent transfer of CH4-derived carbon to viruses. A total of 63% of viral contigs from replicated soil incubations contained homologs of genes present in known methylotrophic bacteria. Genomic sequences of 13C-enriched viruses were represented in over one-third of spacers in CRISPR arrays of multiple closely related Methylocystis populations and revealed differences in their history of viral interaction. Viruses infecting nonmethanotrophic methylotrophs and heterotrophic predatory bacteria were also identified through the analysis of shared homologous genes, demonstrating that carbon is transferred to a diverse range of viruses associated with CH4-fueled microbial food networks.

Metagenomics and Quantitative Stable Isotope Probing Offer Insights into Metabolism of Polycyclic Aromatic Hydrocarbon Degraders in Chronically Polluted Seawater.

Bacterial biodegradation is a significant contributor to remineralization of polycyclic aromatic hydrocarbons (PAHs)-toxic and recalcitrant components of crude oil as well as by-products of partial combustion chronically introduced into seawater via atmospheric deposition. The Deepwater Horizon oil spill demonstrated the speed at which a seed PAH-degrading community maintained by chronic inputs responds to acute pollution. We investigated the diversity and functional potential of a similar seed community in the chronically polluted Port of Los Angeles (POLA), using stable isotope probing with naphthalene, deep-sequenced metagenomes, and carbon incorporation rate measurements at the port and in two sites in the San Pedro Channel. We demonstrate the ability of the community of degraders at the POLA to incorporate carbon from naphthalene, leading to a quick shift in microbial community composition to be dominated by the normally rare Colwellia and Cycloclasticus We show that metagenome-assembled genomes (MAGs) belonged to these naphthalene degraders by matching their 16S-rRNA gene with experimental stable isotope probing data. Surprisingly, we did not find a full PAH degradation pathway in those genomes, even when combining genes from the entire microbial community, leading us to hypothesize that promiscuous dehydrogenases replace canonical naphthalene degradation enzymes in this site. We compared metabolic pathways identified in 29 genomes whose abundance increased in the presence of naphthalene to generate genomic-based recommendations for future optimization of PAH bioremediation at the POLA, e.g., ammonium as opposed to urea, heme or hemoproteins as an iron source, and polar amino acids.IMPORTANCE Oil spills in the marine environment have a devastating effect on marine life and biogeochemical cycles through bioaccumulation of toxic hydrocarbons and oxygen depletion by hydrocarbon-degrading bacteria. Oil-degrading bacteria occur naturally in the ocean, especially where they are supported by chronic inputs of oil or other organic carbon sources, and have a significant role in degradation of oil spills. Polycyclic aromatic hydrocarbons are the most persistent and toxic component of crude oil. Therefore, the bacteria that can break those molecules down are of particular importance. We identified such bacteria at the Port of Los Angeles (POLA), one of the busiest ports worldwide, and characterized their metabolic capabilities. We propose chemical targets based on those analyses to stimulate the activity of these bacteria in case of an oil spill in the Port POLA.

The Functional Significance of Bacterial Predators.

Predation structures food webs, influences energy flow, and alters rates and pathways of nutrient cycling through ecosystems, effects that are well documented for macroscopic predators. In the microbial world, predatory bacteria are common, yet little is known about their rates of growth and roles in energy flows through microbial food webs, in part because these are difficult to quantify. Here, we show that growth and carbon uptake were higher in predatory bacteria compared to nonpredatory bacteria, a finding across 15 sites, synthesizing 82 experiments and over 100,000 taxon-specific measurements of element flow into newly synthesized bacterial DNA. Obligate predatory bacteria grew 36% faster and assimilated carbon at rates 211% higher than nonpredatory bacteria. These differences were less pronounced for facultative predators (6% higher growth rates, 17% higher carbon assimilation rates), though high growth and carbon assimilation rates were observed for some facultative predators, such as members of the genera Lysobacter and Cytophaga, both capable of gliding motility and wolf-pack hunting behavior. Added carbon substrates disproportionately stimulated growth of obligate predators, with responses 63% higher than those of nonpredators for the Bdellovibrionales and 81% higher for the Vampirovibrionales, whereas responses of facultative predators to substrate addition were no different from those of nonpredators. This finding supports the ecological theory that higher productivity increases predator control of lower trophic levels. These findings also indicate that the functional significance of bacterial predators increases with energy flow and that predatory bacteria influence element flow through microbial food webs.IMPORTANCE The word "predator" may conjure images of leopards killing and eating impala on the African savannah or of great white sharks attacking elephant seals off the coast of California. But microorganisms are also predators, including bacteria that kill and eat other bacteria. While predatory bacteria have been found in many environments, it has been challenging to document their importance in nature. This study quantified the growth of predatory and nonpredatory bacteria in soils (and one stream) by tracking isotopically labeled substrates into newly synthesized DNA. Predatory bacteria were more active than nonpredators, and obligate predators, such as Bdellovibrionales and Vampirovibrionales, increased in growth rate in response to added substrates at the base of the food chain, strong evidence of trophic control. This work provides quantitative measures of predator activity and suggests that predatory bacteria-along with protists, nematodes, and phages-are active and important in microbial food webs.

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design.

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets-from soil and seawater microbiomes-to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon's abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.IMPORTANCE One of the biggest challenges in microbial ecology is correlating the identity of microorganisms with the roles they fulfill in natural environmental systems. Studies of microbes in pure culture reveal much about their genomic content and potential functions but may not reflect an organism's activity within its natural community. Culture-independent studies supply a community-wide view of composition and function in the context of community interactions but often fail to link the two. Quantitative stable isotope probing (qSIP) is a method that can link the identity and functional activity of specific microbes within a naturally occurring community. Here, we explore how the resolution of density gradient fractionation affects the error and precision of qSIP results, how they may be improved via additional experimental replication, and discuss cost-benefit balanced scenarios for SIP experimental design.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes.

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Proteorhodopsins dominate the expression of phototrophic mechanisms in seasonal and dynamic marine picoplankton communities.

The most abundant and ubiquitous microbes in the surface ocean use light as an energy source, capturing it via complex chlorophyll-based photosystems or simple retinal-based rhodopsins. Studies in various ocean regimes compared the abundance of these mechanisms, but few investigated their expression. Here we present the first full seasonal study of abundance and expression of light-harvesting mechanisms (proteorhodopsin, PR; aerobic anoxygenic photosynthesis, AAnP; and oxygenic photosynthesis, PSI) from deep-sequenced metagenomes and metatranscriptomes of marine picoplankton (<1 µm) at three coastal stations of the San Pedro Channel in the Pacific Ocean. We show that, regardless of season or sampling location, the most common phototrophic mechanism in metagenomes of this dynamic region was PR (present in 65-104% of the genomes as estimated by single-copy recA), followed by PSI (5-104%) and AAnP (5-32%). Furthermore, the normalized expression (RNA to DNA ratio) of PR genes was higher than that of oxygenic photosynthesis (average ± standard deviation 26.2 ± 8.4 vs. 11 ± 9.7), and the expression of the AAnP marker gene was significantly lower than both mechanisms (0.013 ± 0.02). We demonstrate that PR expression was dominated by the SAR11-cluster year-round, followed by other Alphaproteobacteria, unknown-environmental clusters and Gammaproteobacteria. This highly dynamic system further allowed us to identify a trend for PR spectral tuning, in which blue-absorbing PR genes dominate in areas with low chlorophyll-a concentrations (<0.25 µgL-1). This suggests that PR phototrophy is not an accessory function but instead a central mechanism that can regulate photoheterotrophic population dynamics.

Taxon Disappearance from Microbiome Analysis Reinforces the Value of Mock Communities as a Standard in Every Sequencing Run.

Mock communities have been used in microbiome method development to help estimate biases introduced in PCR amplification and sequencing and to optimize pipeline outputs. Nevertheless, the strong value of routine mock community analysis beyond initial method development is rarely, if ever, considered. Here we report that our routine use of mock communities as internal standards allowed us to discover highly aberrant and strong biases in the relative proportions of multiple taxa in a single Illumina HiSeqPE250 run. In this run, an important archaeal taxon virtually disappeared from all samples, and other mock community taxa showed >2-fold high or low abundance, whereas a rerun of those identical amplicons (from the same reaction tubes) on a different date yielded "normal" results. Although obvious from the strange mock community results, we could have easily missed the problem had we not used the mock communities because of natural variation of microbiomes at our site. The "normal" results were validated over four MiSeqPE300 runs and three HiSeqPE250 runs, and run-to-run variation was usually low. While validating these "normal" results, we also discovered that some mock microbial taxa had relatively modest, but consistent, differences between sequencing platforms. We strongly advise the use of mock communities in every sequencing run to distinguish potentially serious aberrations from natural variations. The mock communities should have more than just a few members and ideally at least partly represent the samples being analyzed to detect problems that show up only in some taxa and also to help validate clustering. IMPORTANCE Despite the routine use of standards and blanks in virtually all chemical or physical assays and most biological studies (a kind of "control"), microbiome analysis has traditionally lacked such standards. Here we show that unexpected problems of unknown origin can occur in such sequencing runs and yield completely incorrect results that would not necessarily be detected without the use of standards. Assuming that the microbiome sequencing analysis works properly every time risks serious errors that can be detected by the use of mock communities.